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ABSTRACT 

We use deep, high quality colour magnitude diagrams obtained with the Hubble Space Telescope to compute a simplified version 
of the Mironov index (SMI; ^) to parametrize the horizontal branch (HB) morphology for 23 globular clusters in the M3I galaxy 
(Sample A), all located in the outer halo at projected distances between 10 kpc and 100 kpc. This allows us to compare them with 
their Galactic counterparts, for which we estimated the SMI exactly in the same way, in the SMI vs. [Fe/H] plane. We find that the 
majority of the considered M3 1 clusters lie in a significantly different locus, in this plane, with respect to Galactic clusters lying at any 
distance from the center of the Milky Way. In particular they have redder HB morphologies at a given metallicity, or, in other words, 
clusters with the same SMI value are » 0.4 dex more metal rich in the Milky Way than in M3 1 . We discuss the possible origin of this 
difference and we conclude that the most likely explanation is that many globular clusters in the outer halo of M3 1 formed « 1 - 2 Gyr 
later than their counterparts in the outer halo of the Milky Way, while differences in the cluster-to-cluster distribution of He abundance 
of individual stars may also play a role. The analysis of another sample of 25 bright M31 clusters (eighteen of them with My < -9.0, 
Sample B), whose SMI estimates are much more uncertain as they are computed on shallow colour magnitude diagrams, suggests that 
extended blue HB tails can be relatively frequent among the most massive M3 1 globular clusters, possibly hinting at the presence of 
multiple populations. 
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1. Introduction 

Globular clusters (GCs) are well studied and widely used tracers 
of the earliest epoch of galaxy formation. They are found in any 
kind of galaxy, they are bright and (typically) several Gyr old, 
at least in the local Universe. Probably the most useful property 
of GCs is that valuable estimates of their fundamental physical 
characteristics, notably metallicity and age, can be obtained with 
various techniques (see Brodie & Strader 2006[ for a recent re- 
view and references). Integrated colours and spectra can be used 
to study unresolved clusters in distant galaxies. In the range of 
distances in which GCs can be resolved into individual stars, 
colour magnitude diagrams (CMD) can be derived and the most 
reliable and precise estimates of the cluster parameters can be 
obtained from observables defined on the CMD. 

This is especially true for cluster age, as the luminosity of 
the main sequence turn off (MSTO) point is unrivaled as an age 
indicator ( [Renzini & Fusi Pecci|1988)|Gallart et al.|2005 ]). The 
advent of large telescopes and modern CCD cameras brought the 
MSTO of all Galactic GCs within reach in the last two decades, 
thus allowing the establishment of a robust and accurate age 



Send offprint requests to: S. Perina 

* Based on observations made with the NASA/ESA Hubble Space 
Telescope, obtained from the Hubble Legacy Archive, which is a collab- 
oration between the Space Telescope Science Institute (STScI/NASA), 
the Space Telescope European Coordinating Facility (ST-ECF/ESA) 
and the Canadian Astronomy Data Centre (CADC/NRC/CSA). STScI 
is operated by the Association of Universities for Research in 
Astronomy, Inc., under NASA contract NAS 5-26555. 



scale (see Dotter et al.||2010 for a state-of-the-art analysis and 
references). Before this epoch, the key age indicator which was 
adopted to obtain global constraints on the Galactic halo forma- 
tion from GCs was the horizontal branch (HB, |Searle & Zinn| 
1978, Zinn 1980). Low-mass stars (m ^ 1 Mq) populate the 
HB in the evolutionary phase of core He burning. The actual 
temperature and luminosity of a given HB star at the begin- 
ning of this phase (Zero Age horizontal branch, ZAHB) is de- 
termined by the complex interaction of several different factors 
(Rood 1973). The overall metallicity is generally recognized as 
the most important parameter, but age. He abundance, stellar ro- 
tation, and any parameter affecting the mass loss along the red 



giant branch (RGB), may have a substantial role ( |Fusi Pecci et 
al 



1993[ [Catelan 2009). The realization that metallicity alone 
was not sufficient to account for the complex behaviour of HB 
morphologies in Galactic GCs led to the so called second param- 
eter problem, that dominated the scientific debate in the field of 



Galactic astronomy fo r a couple of decades (see, e.g., [Sandage 
I& Wildey 1967 ; Zinn||1980| |Preston et al.||1991l |Fusi Pecci et 
al. 1993; Lee et al. 1994; Catelan & de Freitas Pacheco 1994] 



1995,. Fusi Pecci & Bellazzini,|1997> ,Recio-Blanco et al.,,200^ 



Dotter et al.|2 010; Grat ton et al.||2010| and references therein). 



The most recent and thorough analyses concluded that no less 
than three parameters are needed to account for the HB mor- 
phology of Galactic GCs (an idea already proposed in the past, 
see |Fusi Pecci et aL] 1993 Recio-Blanco et al. 2006|l, namely 
metalhcity, age, and He abundance ( Gratton et a .|2010] l. 

The role of He abundance should be considered within the 
framework of the co-existence of multiple populations in glob- 
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ulars (see |Gratton, Carretta, Bragagha|201^ for a recent review 
and references). In the most generally accepted view of the early 
evolution of GCs, subsequent generations of stars have under- 
gone processing of the the light elements, indicative of hot H- 
buming via the CNO cycle, which includes the enrichment of 
He (see, e.g., Carretta et al. 2010 for discussion and references). 



This should lead to strong effects on the HB morphology of 

Caloi &, 



present-day clusters (see, e.g. D'Antona et al. 



^ntona 20q7| |Yoon et al.|2o68||Gratton et al.|2010| ). 



2005, 



While the use of the HB morphology as a possible (albeit 
ambiguos, |Fusi Pecci & Bellazzini|1997) age indicator has been 
superseded by much more reliable techniques for nearby clus- 
ters, it may still be valuable for more distant systems, where the 
MSTO of GCs is too faint to be observed with the instrumenta- 
tion currently available. If we consider the case of our nearest 
neighbour giant galaxy M31, the available observational mate- 
rial on GCs is quite similar to what was available for Galactic 
GCs in the seventies / early eighties, which led to the very influ- 
ential scenario proposed by Searle & Zinn ( 1978| l, 



I.e.: 



- CMDs reaching down to a few magnitudes below the HB 
level are available for several GCs, thanks to Hubble Space 
Telecope (HST) observations (see, e.g. Ajhar et al. 1996 
Fusi Pec ci et al.||19961 [Holland et al.|]T99"71 [Rich et al] 



2005 1 iMackey et al.|2007| [Perina et al.porT Tand references 



therein). It is possible to reach the MSTO luminosity of M3 1 
clusters, but this requires such a large amount of HST time as 
to make a systematic study unpractical ( Brown et al.,2004, ). 
- Global metallicity estimates from integrated spectra and/or 
from the colour of the RGB are also available ( [Galleti et al.| 
[2009; C aldwell et al.|2011||Penna et al.|2009"]l, ana logously 
to the MW GC compilation by [Zinn & West| ( |1984l l. 



In |Federici et al. ( 2012| Paper-I hereafter) we collected a 
dataset of 48 M31 GCs having V,V-I or V,B-V CMD homo- 
geneously derived by us from images collected by different 
group^ In that paper, we used these CMDs to obtain an em- 
pirical calibration of the relation between the absolute V magni- 
tude of the HB (MviHB)) and the metallicity ([Fe/H]), as well 
as a new estimate of the distance to M31. Here we use the same 
dataset, and the cluster parameters homogeneously obtained in 
Paper-I, to get a first systematic outlook on the behaviour of the 
HB morphology in the GC system of that galaxy, in particular in 
comparison with the Milky Way (MW), from the distribution of 
individual stars in each cluster (see Rey et al.|2007 Dalessandro 



|et al .|20l"2| for an independent analysis based on integrated ul 
traviolet colours). We adopt a more straightforward parametriza- 
tion of th e HB m orphologv through a "simplified" version of the 



Mironov 



(1972 1 index MI= 



B+R 



In the following we refer to this 
newly dehned mdex as the Simplified Mironov Index (SMI; see 
Sect. |ZT]for details). 

A previous use of SMI on a smaller sample of M31 GCs 
was made by Rich et al. ( 20051), whereas a more basic approach 



to the HB morphology classification was successfully attempted 
by 'Perina et al. (201 1 1, specifically aimed at verifying the effect 
of the HB morphology on the H^ spectral index. 

The present analysis should be considered as a first attempt 
to provide a view of the M3 1 GC system at a comparable level 
of accuracy as we had about 30 years ago for our own Galaxy 
(but with diff'erent selection biases, see Sect. |2.3|l, which is a 



remarkable step forward in the understanding of the formation 



' The sample includes also four objects classified as Extended 
Clusters (EC), since they appear old and metal-poor as classical globu- 



and early evolution of M31. In Sect.|2]we present our sample and 
assumptions, and we describe the procedure to obtain the HB 
morphology indicator used in the analysis. In Sect.|3]we compare 
the HB morphology of the considered M31 clusters with that 
of MW GCs, and discuss the possible origins of the observed 
differences. Finally, in Sect. |4] we summarize and discuss our 
results also in the context of the literature. 



2. Data and data analysis 

From the dataset presented in Paper-I we selected two samples 
of clusters for our analysis, depending on the quality of the avail- 
able CMD: 



- Sample A: 23 clusters with well defined CMD sequences 
and easily identifiable HB, limiting magnitude Vq ^ 27.0 
(i.e. reaching approximately 2 mag below the HB level, at 
least), and low contamination from the M3 1 field. This con- 
stitutes the bulk of our analysis. The CMD for all the clus- 
ters in Sample A are shown in Fig. [T] to give the reader a 
clear idea of the quality of the HB morphology classification 
that can be obtained from this material. It must be noted that 
the quality of the data for Sample A clusters is not fully ho- 
mogeneous. The CMDs of B384 and B468 have a brighter 
limiting magnitude than most of the other clusters of this 
sample. B008, B292, B298, B336, B337, B350, and B531 
suffer from a slightly higher degree of contamination and/or 
crowding, w.r.t to the other Sample A clusters, and/or may 
be affected by some differential reddening (still not seriously 
affecting their HB morphology). 

- Sample B: 25 clusters with CMD not fulfilling the above cri- 
teria, still valuable, however, to obtain an estimate of SMI, 
albeit significantly more uncertain than for Sample A. We 
briefly consider Sample B here only for completeness and to 
get preliminary insight for future investigations. 



In Paper-I all the cluster CMDs were compared with a grid of 
RGB and HB templates of Galactic GCs with a fitting procedure 
that provides simultaneous best estimates of the HB magnitude 
level Vhb, the reddening E(B-V), the metallicity [Fe/H], and the 
distance modulus for each cluster (see Paper-I for a detailed de- 
scription and discussion about uncertainties; we discuss the role 
of these uncertainties in our analysis in Sect.|3]l. In the following 
we adopt Vhb, E(B-V) and [Fe/H] from Paper-I; it is clear that 
Sample A clusters, having the best CMDs, have also the most 
reliable and accurate estimates of these parameters. 

All the results of the present study are based on the compari- 
son with Galactic GCs taken as a reference. As an homogeneous 
Galactic sample, we adopt the publicly ava i lable HST photom- 



etry of 74 Galactic GCs from 
of [Fe/H], Vhb, E(B-V) and (m 



Piotto et al. 



( |2002|p 



The values 
M)o for these clusters have 
been taken from Harris ( 1996 2010 version). The adopted ex- 
tinction laws are the same as in Paper-I, Ay = 3.1E{B - V), 
Ai = 1.94E(B - V) and E(V - /) = 1.375£(B - V). The Ust of 
Sample A and Sample B clusters and the relevant parameters for 
the present analysis are provided in Table [T] 



lars (see Mackey et al. 20061 and references therein) 



- The original data jPiotto et al.|2002| were obtained in the F439W 
and F555W bands of the Wide Field and Planetary Camera 2 (WFPC2). 
Photometry in these filters has been converted into standard B,V pho- 
tometry as described in Paper-I. 
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Fig. 1. colour-magnitude diagrams of the target clusters. The colour on the x axis is indicated by a label in the top-right corner of 
each panel : BV corresponds to the dereddened colour (B - V)o and VI corresponds to the dereddened colour {V - I)o- 



2.1. The Simplified Mironov Index (SMI) 
In the widely used original Mironov Index MI=^|^(see 



e.g. [Searle & Zinn|1978HLee et al.|1994}|Catelan & de Freitas 
|Pacheco| 1994 1 and references therein), R and B are the numbers 
of HB stars respectively redder and bluer than the RR Lyrae in- 
stability strip edges. Here we use the same formulation, but R 
and B are defined as the number of HB stars redder and bluer 
than a given colour threshold, approximately located in the mid- 
dle of the RR Lyrae instability strip. Given the sometimes scanty 
population on the HB, this appeared as the most sensible choice 
allowing us to account for all the HB stars with a single param- 
eter. 

We estimated the SMI from extinction-corrected CMDs as 
illustrated in Fig.|2]and described below: 

1 . For each cluster we consider a circular area around the center 
that is clearly dominated by cluster members. All the star 
counts described below are performed in this selected circle. 

2. We count as R the HB stars lying within ±0.5 mag from Vhb 
and having 0.50 < iV-I)o < 0.80or0.30 < (B-y)o < 0.65, 
depending on the photometric bands of the available CMD. 

3. We count as B the HB stars having {Vhb - 1.0) < Vq < 26.0, 
and {V - I)o < 0.50 or {B - V)o < 0.3, depending on the 
available CMDs. The Vq < 26.0 limit has been adopted to 
avoid regions of the CMD were the completeness level can 
be sharply falling, thus strongly biasing the star counts. Note 
that the same selection boxes are adopted for M31 and MW 
GC^ This implies that a fraction of genuine and clearly 



The Vo < 26.0 limit for the selection of B stars adopted for M31 
clusters correspond to My = 1.58, adopting the average distance mod- 
ulus from Paper-I, (m - M)o = 24.42. This has been converted into 
the corresponding limit for each Galactic GC adding their true distance 
moduli. 



identifiable Blue HB (BHB) stars can be excluded in some 
Galactic GCs (as in the case of NGC6229, shown in Fig. 2]l, 
to preserve the full homogeneity in the comparison with M3 1 
clusters. The case of NGC6229 shows also that the same 
threshold helps to prevent contamination of the BHB star 
counts from blue stragglers. In any case, stars appearing to 
be associated with the Blue Straggler sequence are excluded 
from the computation of SMI. 

For Sample A clusters, the average number of field stars 
falling in each selection box was estimated using an annulus 
(field annulus) surrounding the cluster with the same area as 
the selected circle. Field counts are subtracted from the clus- 
ter counts and SMI is computed propagating Poisson errors 
on the star counts. The effects of uncertainty in Vhb are neg- 
ligible, because of the generous interval in V that is consid- 
ered in the selection of HB stars. The uncertainty in E(B-V) 
can have a significant impact, in principle: we consider the 
case in detail in Sect. [3] 

For Sample B clusters, we directly estimate B and R from 
statistically decontaminated CMDs (see |Perina et al.||2009] 
and Paper-I). This is a less robust procedure for star counts 
with respect to that adopted for Sample A, but it was the only 
viable one, since the identification of key CMD features was 
too uncertain without previous decontamination. 



The location of the threshold between R and B stars is arbi- 
trary and has been set only for convenience. However, since the 
aim is the comparison between two sets of GCs where the SMI 
has been computed in the same way, this cannot introduce any 
serious bias in the analysis. On the other hand, the consistency 
between the threshold in {B - V)q and {V - /)o has been checked 
for accuracy using clusters that have photometry in all the three 
passbands. 
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Fig. 2. Illustration of the selection boxes adopted for the compu- 
tation of the Simplified Mironov Index. Left panel: CMD of the 
Galactic globular cluster NGC6229, reported to the mean dis- 
tance of M3 1 . Right panel: CMD of the Sample A globular clus- 
ter MCGC5. In each panel, the vertical segments mark, from left 
to right, the boundary between the blue (B) and red (R) part of 
the HB and the red limit of the HB selection. The continuous 
horizontal line is the mean level of the HB, while the dotted hor- 
izontal segments enclose the stars selected as bona-fide HB. In 
both diagrams Blue and Red HB stars selected to compute SMI 
are plotted in dark and light grey, respectively. 



In a few clusters the red HB is totally or partially superposed 
to the RGB. In these cases we estimated the contribution of these 
stars to R by constructing the luminosity function (LF) of RGB 
stars in the selected circle and in the corresponding field annu- 
lus. Once subtracted the latter from the former, to remove any 
possible contamination from field Red Clump stars, we look for 
residual peaks in the decontaminated LF that can correspond to 
the red HB. Great care is taken to discriminate between the red 
HB and the (fairly smaller and, generally, brighter) peak pro- 
duced by the RGB bump. Then the number of stars in the red 
HB peak is estimated by subtracting the underlying RGB, inter- 
polated in the region of the red HB peak. We checked several 
clusters in this way for the presence of red HB stars superposed 
to the RGB, and a significant signal was detected only for two 
clusters, MCGC8 (36 R stai-s) and MCECl (46 R stars). In some 
cases, like e.g. B468, although the red HB is clearly superposed 
to the RGB, the above procedure is pointless as SMI vanishes all 
the way. For B298 and MCEC3 the uncertainties on the number 
of R stars are somewhat larger because of a stronger field con- 
tamination. For these two clusters we provide only upper limits 
to their SMI, but the true value cannot be much different from 
the reported one. 

2.2. Relation with other HB morphology indicators 

It is clear that the SMI should have a similar behaviour as the 
original MI which, in turn, is a simpler version of the more gen- 
erally used horizontal branch Ratio HBR= bIrIv (where V is the 



number of stars in the RR Lyrae instability strip, Lee et al. 1994| l. 
The two indices are tied by the relation ( [Preston et al.|199T] i: 



D D D 

= 0.50 + 0.50- 



B+R 



B + R + V 



(1) 



Fig.|3^ shows that the same relation provides an excellent fit 
also to the trend between SMI and HBR. 

These indices trace essentially the peak (or the weighted 
mean) of the HB colour distribution, this being especially true 
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Fig. 3. Panel (a): Correlation of the Simplified Mironov Index 
with the HBR index (gf^fy) for our sample of Galactic GCs. 
HBR values are taken from the latest version of the Harris ( 1996]l 
catalogue. The continuous Une is the relation between ^^and 



from 'Preston et al. ( 1991 1, reported in Eq. [T] Panel (b): 



Correlation of SMI with trie integrated ultraviolet-optical colour 
{NUV-V)o for Sample A (open stars). Sample B (open triangles, 
to be considered as lower limits), and Galactic GCs (small filled 
circles). 



for the SMI because the Vq < 26.0 limit makes it blind to stars 
in the extreme blue tail of the HB. Other parameters - e.g. L, or 
T"'?f(HB) - are better suited to trace the extension of the HB, in 



eff ,——11 

particular of the blue tail, when present (Preston et al . ' 1 99 1 'Fusil 
Pecci et al. 1993; Recio-Blanco et al. 2006; Gratton et al. 20101^ 



It has been shown ( |Rey et al. 2007| [Dalessandro et al.| 
|2012j ) that the integrated near/far ultraviolet - visible colours 
(NUV - V)o and (FUV - V)o are good tracers of the GC HB 
morphology. In principle they should be preferred to our SMI, 
because by definition they cannot be affected by the incomplete- 
ness problem that forced us to adopt the Vq < 26.0 limit, and 
so are particularly sensitive (especially (FUV - V)o) to the blue 
tails of the HB distributions (Dal essandro et al.|20I2| l. However, 
these integrated colours are much more affected by reddening 
and contamination from hot non-HB stars (e.g. from the M31 
disk field population or the turn-off stars of the GC itself, espe- 
cially {NUV-V)o; see, e.g., Sect. [ZZT] ). Therefore, SMI should 
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provide a complementary view of the HB morphology, with re- 
spect to the UV-optical integrated colours. 

In Fig.jSj) we show that indeed SMI and {NUV -V\) broadly 
anti-correlate, and the behaviour of Galactic and M31 clusters 
from our samples are very similar in this plan^ We note that 
most of the M31 clusters in our samples having available NUV 
and FUV magnitudes are from Sample B. It is reassuring to note 
that, in spite of the large uncertainties associated to the SMI esti- 
mates for Sample B clusters, SMI anti-correlates quite well with 
iNUV-V)Q. The same degree of correlation is observed between 
(NUV - V\) and HBR (see e.g. Dalessandro et al. 2012). 

Given the results illustrated in Fig. |3] we consider that our 
SMI is validated as a reliable HB morphology indicator tracing 
the peak of the HB distribution as MI and HBR. For most of the 
clusters in Sample A which lack UV-optical colours, this is the 
only parametrization of the HB morphology currently available. 

2.2.1 . SMI vs. (FUV - y)o: the effect of blue tails 

The comparison between SMI and (FUV -V)(i, shown in Fig. 4, 
deserves a deeper discussion. The anti-correlation is signifi- 
cantly less tight than in the SMI vs (NUV -V\) plane and, above 
all, the distribution of M3 1 clusters is different from their MW 
counterparts. In particular they are all confined in a relatively 
narrow strip along the blue edge of the (FUV — V)q distribution. 

This is partially due to the selection in FUV flux that affects 
the M31 sample. At that distance clusters with FUV fluxes be- 
low a certain threshold were not detected by GALEX: in practice 
M31 GCs with Mv ^ -7.0 do not have valid measures of the in- 



tegrated FUV magnitude ( Rey et al.|2007| l. This translates into a 
threshold in (FUV - y)o colour ranging from (FUV - V\) $7.2 
for the reddest clusters (having integrated (B - y)o ~ 1.0) 
to (FUV - y)o 5 5.0 for the bluest ones (having integrated 
(B — y)o ~ 0.3), effectively squeezing the M31 GCs in our sam- 
ple into a narrower range of (FUV — y)o with respect to MW 
clusters. 

However, there are two additional effects that can concur 
to produce the observed difference. First, SMI estimates for 
Sample B clusters (the large majority of M31 clusters in Fig. 4) 
should be considered as lower limits, because the bright limiting 
magnitude of their CMD may wipe out most (or all) of their 
BHB stars (if present) from the computation of SMI. Larger 
SMI values for some of these clusters would help to distribute 
them over the same range covered by Galactic GCs in this plane. 
Furthermore, some of them may have very extended blue tails, 
including a population of extreme HB (EHB) and/or (possibly) 
blue hook stars (see |Catelan||2009| for references and discus- 
sion). 

The inset CMD for three Galactic clusters (NGC6584, 
NGC1851 and NGC2808) having similar SMI and widely dif- 
ferent (FUV - y)o illustrates this possibility and shows very 
clearly the complementarity of SMI and (FUV - V\) in describ- 
ing the HB morphology. All the three clusters have a bimodal 
HB morphology, however they greatly differ in the extension of 



* Integrated FUV and NUV magnitudes (from Galaxy Evolution 
Explorer - GALEX photometry) are taken from Re y et af] ( |2007[ l, for 
M31 clusters, and from Dalessandro et al. (2012) for MW clusters. The 
values of (NUV - V)o and (FUV - V)o have been computed adopting 
AA,t,v=8.90E(B-V) and Aft,v=8.16 (from Rey et al. 2007) , and taking 
V magnitudes from RBC V4.0 and E(B-V) from Paper-I. The adoption 
of the very recent updated and extended dataset by |Kang et aI.| ( [20T2) 
would not lead to any significant improvement, for the sample of M3 1 
clusters considered here. 



their blue tail. The hottest temperature reached by their HB cor- 
responds to the largest integrated FUV flux and, consequently, 
the bluest (FUV - V)q colour (see |Dalessandro et al.|2012 ). The 
difference in the morphology between NGC4590 and NGC6229 
(also labelled in the figure and discussed in Sect. |3] below) is 
similar to that between NGC6584 and NGC1851/NGC2808. 

The trend shown in the CMDs indicates that a remarkably 
extended BHB is required to reach the bluest limits of the 
(FUV -V)(i range, thus suggesting that several sample B clusters 
likely have a significant blue tail population that went undetected 
in their shallow CMDs. Hence extreme HB morphologies like 
that displayed by NGC2808, and a few other peculiar clusters in 
the Milky Way, i.e. to Cen, M54 and NGC2419, may be relatively 
common among Sample B clusters (see also [Dalessandro et aL] 
!2012). Interestingly enough, such extreme (FUV - V)o colors 
indicate that a substantial fraction (3-6%, |Buzzoni et al.|[20T2) l 
of the total (bolometric) luminosity for these clusters is emitted 
shortward of 2500 A, a feature that may closely deal, on a larger 
galactic scale, with the striking phenomenon of the UV upturn, 
as extensively observed among elliptical galaxies ( |Yi & Yoon| 
|2004l ). 

The case of B337 is especially interesting since it has 
SMI=0.0, i.e. purely red HB as seen from its CMD, and ex- 
tremely blue (FUV - V)o colour Still, it is a sample A cluster 
and in Fig. [T|it can be appreciated that no significant BHB pop- 
ulation emerges in the CMD down to » 3 mag below the red HB 
level, a range containing the majority of BHB stars in NGC2808. 
A detailed inspection of the available HST images revealed the 
presence of two possible contaminating sources lying within 
the aperture adopted by |Rey et al.| ( |2007| l for their NUV/FUV 
GALEX photometry (r„p = 4.5")- A star nearly two magnitudes 
brighter than the brightest cluster star is located at x 3.5" from 
the cluster center Its colour ((B-V)o - 0.45) and magnitude are 
typical of a foreground Galactic dwarf. While this bright star is 
bluer than most of the cluster giants and clearly contaminates the 
integrated FUV magnitude of the cluster, it is too red to push the 
integrated (FUV - V)o to the extreme value observed. Another 
potential contaminant has been identified at ^ 0.25" from the 
cluster center, a region where photometry packages are unable 
to reliably resolve and estimate fluxes of any star. A very rough 
estimate based on the comparison of the intensity peak with stars 
of known magnitudes suggest that this star is brighter than the 
cluster RGB tip and may have a color as blue as (B - V)o « 0.0. 
This would be consistent with a Post Asymptotic Giant Branch 
star (Jasniewicz & Parthasarathy 2009), that may well dominate 
the UV flux of the cluster. Hence, B337 may be the example 
of a cluster whose integrated UV-optical colour is strongly af- 
fected by fore/back-ground contamination and/or by a rare non- 
HB source. However, these cases should be quite unusual in our 
samples. 



2.3. An obvious selection bias 

Sample A is clearly not representative of the whole M3 1 GC sys- 
tem. The quality criteria imposed on the CMD are very hard to 
fulfill for clusters that are projected on the high surface bright- 
ness bulge or disk of M31, since there (a) the background level 
and the crowding are high, preventing accurate photometry of 
individual stars down to the required faint magnitude limit, and 
(b) the contamination by field M3 1 stars can be also very high, 
strongly affecting star counts on the HB. Moreover,clusters more 
distant from the center of their parent galaxy have intrinsically 
larger sizes ( |van den Bergh|[T994) ), hence remote M31 clusters 
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Fig. 4. Simplified Mironov Index versus the integrated ultraviolet-optical colours (FUV - V)o for Sample A (open stars). Sample B 
(open triangles, to be considered as lower limits) and Galactic GC (small filled circles). Three remarkable Galactic GCs are high- 
lighted with a concentric open circle, are labelled in boldface and have their CMD presented in an inset box, to show the actual 
variation of their HB morphology with decreasing {FUV-V)o at nearly constant SMI. Two other Galactic clusters and one Sample A 
cluster, that are mentioned in the text, are also labelled. 



are less severely crowded than those near the center, in average, 
when seen from the Earth. 

Indeed, the left panel of Fig.|5]shows that most of Sample A 
clusters lies at projected distances Rp > 10 kpc from the cen- 
ter of M31, a radius that encloses approximately 80 per cent of 
the whole sample of M3 1 GC^ Moreover, ^ 20 per cent of 
Sample A clusters have Rj, > 40 kpc, a range that contains less 
than 5 per cent of the whole sample. Finally, many of the outer- 
most clusters of Sample A have been found to be likely associ- 
ated with known existing substructures in the outer halo of M3 1 
( [Perina et al^[2009l [Mackey et aT1[20T0| ). Therefore, Sample A 



^ The data are from the Revised Bologna Catalogue of M31 
globular clusters and candidates [Galleti et aT] ( |2004[ RBC V4.0); 
http : /www . bo . astro . it/M3 1/ 



is more likely representative of the outer halo population of the 
GC system of M31, and/or of the population of GCs more re- 
cently accreted via the disruption of their parent dwarf galaxies. 
This selection bias should be kept in mind very clearly when 
comparing Sample A with Galactic globulars. 

Sample B shows a much more concentrated distribution of 
Rp, more similar to the overall population. However, the right 
panel of Fig. |5] shows that the distribution of spatial (3-D) dis- 
tance^is very similar to the one of Sample A, and clearly also 
not representative of the overall M3 1 GC system. 



* Derived by combining the projected distances and the line-of-sight 
distances from the center of M3 1 , the latter obtained from the individual 
distance moduli from Paper-I, as done also by [Mackey et al.| j2006| 
2007). 
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Fig. 5. Left panel: cumulative distributions of projected dis- 
tances from the center of M31 of Sample A (long dashed line) 
and Sample B (dotted line) are compared with the distribution 
of all the confirmed M3 1 clusters classified as old (GCs) (solid 
line). Right panel: the same comparison but for the distributions 
of three-dimensional distances from the center of M31 (Rcc)- 

For the whole sample we used J^Rp as a proxy for Rcc- 



recent 



It is worth recalling that Mackey et al. ( 2006[ l and Mackey 



2007| l noted a strong second parameter effect among the 



y discovered group of remote (and/or extended) clusters 
of M 31 (Hu xor et aL 12004] |2005l [GaUeti e t al. 2006 ). However, 
a quantitative and systematic comparison with MW clusters is 
attempted here for the first time. 



3. Results 

The main result of the present analysis is shown in Fig. |6] 
Sample A clusters are compared with MW clusters in the clas- 
sical plane opposing metallicity and HB morphology. The solid 
lines show three theoretical isochrones illustrating the expected 
trend of the SMI as a function of [Fe/H]. These have been ob- 
tained by converting isochrones in the ;gf^fy vs. [Fe/H] plane 
(derived by Rey et al. 2001 from synthetic HB populations mod- 
els) using Eq.yj 

It must be stressed that these isochrones are plotted only as 
a general reference for the interpretation of the observed distri- 
butions, since they (a) refer to a slightly different quantity than 
that plotted in Fig.|6](i.e. the original MI instead of SMI, but this 
does not seem a reason of serious concern, see Fig. [3]l, and (b) 
the details of their shape depend on the mass loss recipe adopted 
in modeling the synthetic H B (Rey et al. 2001; |Lee et al.|1994[ 
[Catelan & de Freitas Pacheco.1994) . The three isochrones cor- 
respond to different ages in steps of 1. 1 Gyr, setting a rough rel- 
ative age scale, from top to bottom At=0, -1.1 and -2.2 Gyr (see 
Fig. 9 in Rey et al. 2001). Alternatively, they can be read as a 



tionQ We stress that, for the same reasons mentioned above for 
the age scale, this should be considered just as a rough indica- 
tion of the effect Y on SMI. It is also important to keep in mind 
that in our current view of globular clusters, a previously usual 
concept like "the He abundance of a given cluster" is no more 
valid, since all GCs presenting a spread in light elements are ex- 
pected to present a spread also in Y (see Sect. [T]). Hence each 
of them is characterized, in this respect, by the distribution of 
He abundances of its member stars. Different distributions of Y 
would lead to differences in HB morphology and, consequently, 
in SMI. However both the average abundance and the size of the 
spread (as well as the shape of the distribution) would concur 
in determining the final SMI value of a given cluster. Hence any 
difference in SMI caused by helium, should trace a difference in 
the Y distributions of the considered clusters. 

Fig. [6 shows the well known fact that clusters residing in 
the outer halo (OH, Rcc > 8.0 kpc) of the MW tend to lie on 
an isochrone corresponding to a younger age with respect to 
clusters in the Inner Halo (IH, Rcc < 8.0 kpc |Lee et al.|1994 



Catelan & de Freitas Pacheco|1994|[T995||Rey et al.|200lf . This 
is especially evident at intermediate metalUcity, since at the ex- 
tremes of the metallicity range the HB morphology parameters 
like SMI saturate (i.e. are not able to discriminate within HB 
distributions that have only red or only blue stars). The recent 
analyses by|Gratton et al.| (|2010| l, [Potter et ak] ( |2010| l, and, in 
particular, Dotter et al. (2011 1 confirm that the observed differ- 
ence largely traces actual age differences between the two sub- 
groups of Galactic GCs: while IH clusters are approximately co- 
eval at any metallicity, several OH clusters with [Fe/H]^ -1.5 
are ^ \ - 2 Gyr younger than their IH counterparts with the 
same metallicity. 

The completely new feature of Fig. |6]is the location of M31 
GCs: most Sample A clusters with [Fe/H]<-1.2 lie in a different 
locus with respect to both IH and OH Galactic GCs, approx- 
imately along the isochrone corresponding to At — —2.2 Gyr. 
Specifically, the clusters lying near (or below) the Af - -2.2 Gyr 
isochrone are: B337, MCGC9, MCGC8, MCEC2, MCGC5, 
MCEC4, B336, B514, MCGCIO, MCECl, MCGCl. Therefore, 
eleven out of seventeen Sample A clusters in the considered 
metallicity range, i.e. more than half the sample, lies near this 
locus. These GCs (all labelled in Fig. |6]l appear to have sig- 
nificantly redder HB morphology with respect to their Galactic 
counterparts of the same metallicity. From a different perspec- 
tive, one can say that a given SMI value is reached by these M31 
clusters at a lower metallicity by « -0.4 dex than their Galactic 
counterparts. In the following we will refer to these eleven M31 
clusters as Anomalously Red (HB) Sample A clusters, abbrevi- 
ated with the acronym ARSA. Note also that, for [Fe/H]:$ -0.8, 
no Sample A cluster lie near the At — 0.0 isochrone, along which 
most of the Galactic GCs are clustered. 

Only two Galactic GCs are lying on the Af = -2.2 Gyr 
isochrone as the ARSA clusters, i.e. NGC 4590 and NGC 7078. 
Both clusters are very metal-poor ([Fe/H]:$ -2.0) and are not far 
from the blue saturation limit of the SMI. The inspection of the 
CMD of NGC 7078 reveals that the faint limit imposed by our 
definition of SMI (V < 26.0 for M31 clusters) cuts out a large 
number of blue HB stars belonging to the extended blue tail. 



relative scale of He abundance Y, approximately corresponding 
to AY=H-0.019, 0.0 and -0.019 (M. Catelan, private communica- 



' The adopted standard Y abundance is Y=0.23-l-2Z, where Z is the 
fractional abundance of heavy elements l |Yi et al.|2001^ . See Fig. 9 of 
|Rey et al"| (2007 ) for an example of how the shape of this kind of HB- 
morphology isochrones can be changed by the assumption of a signifi- 
cant spread in He abundance within the HB population. In our case, the 
presence of a population of He-rich stars would prevent the isochrones 
from reaching SMI values near zero even at high metallicities. 
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namely ^30% of the whole HB population. Therefore, in this 
case the SMI does not provide a good parametrization of the HB 
morphology, which is actually bluer than indicated by the SMI. 
We note, however, that virtually all the blue HB stars excluded 
from the SMI estimate for NGC 7078 he within ^ 1.5 mag of 
the threshold (similar to the case of NGC 6229, illustrated in 
Fig. |2]l. Fig. [T] shows that none of the ARSA clusters seems to 
have such a large population of blue HB stars in this range of 
magnitude. So it is very unlikely that the SMI of these clusters 
can be severely affected by this kind of problem. 

On the other hand, no HB star of NGC 4590 lies below 
the faint limit of the selection, and hence the SMI accounts for 
the whole HB distribution, and its relatively low value traces a 
real difference in HB morphology between this cluster and other 
GCs of similar metallicity. This peculiarity is remarked also by 
[Dalessandro et al.| ( |2012| l, who showed that there is a deficiency 
of very blue HB stars in this cluster with respect to, for exam- 
ple, NGC 7099. Since NGC 459 appears to be as old as the 
other metal-poor GCs of the MW, [Dales sandro et al.| ( p0r2,) sug- 
gest that the observed difference in HB morphology may be ac- 
counted for by a difference in He abundanca^ However, it has to 
be recalled that the apparent anomalous position of these clus- 
ters may also be due to inadequacies of the theoretical models, 
especially in this low metallicity regime, where the shape of the 
isocrones may strongly depend on the assumptions on mass-loss 
(see Zinn|1993 Catelan|2009 for examples and discussion). 



Obviously, the unusual location of an individual M3 1 GC in 
the SMI vs. [Fe/H] plane may be considered as hardly significant 
given the large uncertainties. However, to appreciate the actual 
sensitivity of SMI for the best Sample A clusters it may be inter- 
esting to consider the cases of MCGC5 and B514. Both clusters 
have very deep and clean CMDs, they have the same metallicity 
and they differ by 0.45 + 0.13 in SMI, formally a a: 3.5cr differ- 
ence. The comparison of their CMDs (in Fig.[T]i reveals at a first 
glance that the HB morphology of B5 14 is indeed much bluer 
than that of MCGC5, even if both clusters have stars to the red 
and to the blue of the SMI B-to-R threshold. The same visual 
comparison can be made for all the ARSA clusters: the cluster- 
to-cluster differences measured by SMI corresponds to real dif- 
ferences in HB morphology that, in most cases, can be appreci- 
ated simply by looking at the CMDs. Moreover, Sample A clus- 
ters clearly display a collective behaviour that is different from 
MW GCs (and similar to dwarf galaxy satellites of the MW and 
M3 1, see Yang & Sarajediru 2012 1, thus strongly hinting at a real 
difference in some fundamental physical parameter between the 
samples. 



3. 1 . The effect of systematics 

Since the observed difference seems to involve the large majority 
of Sample A clusters, it would be easily explained by a system- 
atic error in one of the relevant parameters, i.e metallicity and 
reddening, or SMI itself. 



* We note that the two couples of clusters considered by'Dalessandro] 
[etai., (2012), i.e. NGC4590 vs. NGC7099 and NGC5466 vs. NGC6341, 
display also a large difference in My and central density which may 
play a role in the chemical enrichment of cluster stars. The correlation 
between these parameters and the extension of blue tails has been noted 
and discussed in Fusi Pecci et al. i 1993 1 and Recio-Blanco et al. j2006| . 



3.1.1. SMI 

The most obvious candidate for a systematic effect is our HB 
morphology parameter If large numbers of BHB stars were 
missed in our computation, the actual SMI would be mere lower 
limits, as in the case of Sample B (see Sect. 3.2 1. Once taken into 
account, these hypothetical BHB stars would lead ARSA clus- 
ters to move near the MW clusters, at least those residing in the 
outer halo of the MW. 

However, the inspection of the CMD of ARSA clusters 
(Fig. TJ reveals that there is no sign of such a significant addi- 
tional BHB population down to a; 3 magnitudes below the Vhb 
level. Thus, any hidden BHB should be made by extreme HB 
and/or blue hook stars, like the faintest HB stars in NGC2808, 
i.e. stars that would be excluded from the computation of SMI by 
definition and cannot be at the origin of the difference between 
M31 and MW GCs shown in Fig.|6] Moreover, the presence of a 
significant population of such stars after a gap in the HB distri- 
bution of ^ 2 mag would make the overall morphology of ARSA 
clusters exceedingly peculiar, calling for an interpretation as the 
results discussed here. 



3.1.2. Metallicity 

A systematic underestimate of the metallicity of M31 clus- 
ters by X 0.4 dex would move all the ARSA clusters on the 
same isochrone populated by Galactic OH clusters. However, 
such a large systematic error can be safely excluded. First, the 
metallicity estimates used here have been obtained in Paper-I 
from the comparison of the observed RGB with templates of 
Galactic GCs, and hence the adopted metallicity scale is ex- 
actly the same for the MW and the M31 samples (except for 
a possible factor that is discussed below). Second, in Paper-I 
we have compared our [Fe/H] estimates with three sets of in- 
dependent spectroscopic estimates and we found that system- 
atic differences are at most 0.1-0.2 dex and in any case con- 
sistent with zero (see Fig. 6 of Paper-I). It is particularly reas- 
suring that the agreement with estimates from Fel lines from 
high resolution spectroscopy by Colucci et al. ( 2009 | l is excel- 
lent iA[Fe/H]cMD-F<.i = -0.04 ± 0.21). 

An implicit assumption in the metallicity scale of Paper-I 
is that M31 clusters are as enhanced in a-element abundance 
(w.r.t the Sun) as their Galactic counterparts. This assumption 
was found to be valid, at least for the handful of bright M31 
clusters analyzed by jColuccTet al.]pO09| ). However this is not 
necessarily true for Sample A clusters. If these clusters had a 
solar [a/Fe] ratio, the comparison with the grid of templates of 
Qf-enhanced Galactic GCs would lead to underestimate the actual 
metallicity. Using the formulae provided by Ferraro et al. ( I999| l 
we estimated that in this case the [Fe/H] values of Sample A 
clusters should be increased by -hO.2 dex to keep them in the 
same metallicity scale as the Galactic GCs. This would reduce 
the difference between the M31 and the MW clusters shown in 
Fig. [6] but it would not be sufficient to cancel it out, since the 
observed effect is 2 times larger Even if a difference in [a/Fe] 
were responsible for part of the observed effect, this would be 
nevertheless an important difference in a fundamental physical 
parameter of the clusters (i.e. the [a/Fe] ratio), revealed thanks 
to Fig. [6] 

3.1.3. Reddening 

The other parameter that can produce a systematic effect is red- 
dening. Like metallicity and distance, reddening was estimated 
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Fig. 6. The SMI vs. metallicity diagram. The Galactic GCs are plotted as filled circles, the M31 Sample A GCs as open stars. 
The symbol colour code shows the galacto-centric distance of the cluster, as described in the bottom-left corner of the panel. The 
Rcc - 20 kpc threshold was somehow arbitrarily adopted for Sample A clusters to highlight those that cannot be considered as 
exceptionally far, since it correspond to the radius containing » 90% of the whole population of confirmed M31 GCs (see Fig. |5]l. 



Arrows indicate upper limits. Solid lines are isochrones from the synthetic HB models by Rey et al. (2001 1, and are labelled 



according to their age difference in Gyr Typical error-bars for M31 GCs are shown in the top-right corner Sample A clusters lying 
near the At - -2.2 Gyr isochrone are labelled. 



simultaneously in Paper-I by a best-fit of the CMD with the grid 
of MW GC templates, and it has been verified that the best-fit 
value cannot be changed by more than +0.02 mag without com- 
promising the overall fit by obtaining unrealistic (or nonsense) 
values of distance and metallicity. Moreover, changing the red- 
dening has little direct effect on the final SMI value (by shifting 
the Blue/Red threshold), while it may have a sizable effect on 
the estimated metallicity. It turns out that to obtain the higher 
[Fe/H], which would be necessary to reduce the difference be- 
tween M31 and MW clusters in Fig. [6] a lower reddening must 
be adopted with respect to Paper-I. However, the comparison of 
our reddening values with three independent sets of estimates 
demonstrates that, if any small systematic effect were there, it 



would go in the opposite sense, since our E(B-V) are on average 
a few hundredths of a magnitude smaller than those estimated 
by other studies (see Paper-I for details). 

However, for the sake of completeness, we have repeated, for 
the ARSA clusters, the whole procedure of metallicity estimate 
performed in Paper-I and the computation of SMI, but forcing 
lower reddening values by 0.02 mag with respect to the best-fit 
solutions of Paper-I. It turns out that (a) the changes in SMI are 
negligible (< 0.05), and (b) [Fe/H] values typically decrease by 
0. 1 dex, with a range between 0.05 dex and 0.25 dex, in agree- 
ment with the predictions by Buzzoni (1995] his Eq. 19). While 
this clearly helps in reducing the difference between Sample A 
and the MW sample, it is not sufficient to close the gap, and for 
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Fig. 7. Same as Fig.|6]but for Sample B M31 clusters. All the es- 
timates should be considered as lower limits, since the bright 
limiting magnitude of the available CMD does not allow to 
reveal a Blue Tail HB component, if present. The metal-rich 
Galactic GCs displaying bimodal HB morphology NGC6441 
and NGC6388 are indicated with labels. Sample B clusters dis- 
cussed in the text are also labelled. 



this 'partial result' we pay a toll: all the CMD fits obtained with 
the lower E(B-V) values are significantly worse than what found 
in Paper-I, and in most cases clearly not acceptable. 

In conclusion, it seems very unlikely that a systematic mises- 
timate of [Fe/H] or E(B-V) can be the (only) cause of the differ- 
ent behaviour of Sample A and MW clusters reported in Fig|6] 

3.2. Sample B 

In Fig. [T] we show the distribution of Sample B clusters in the 
SMI vs. [Fe/H] plane. Unfortunately, the uncertainties are too 
large to draw any firm conclusion, as most of the clusters in the 
interesting range of metallicity have just rough estimates of their 
SMI. In particular, all the SMI estimates for Sample B clusters 
should be considered as lower limits, since the bright limiting 
magnitude of the available CMDs does not allow to reveal even 
the reddest portion of the BHB component, if present. However, 
the comparison with other independent tracers of the HB mor- 
phology shown in Fig. 4 demonstrates that the estimated SMI 
may still carry useful information on the HB of these clusters. 

Taking Fig. |7] at face value, it can be noted that several 
Sample B clusters lie on the isochrones that are populated also 
by MW GCs, even at intermediate metallicity. Still, there are a 
few potentially interesting outliers: 

- four anomalous clusters lying on the Af - -2.2 Gyr 
isochrone, as most Sample A clusters do, namely B027, 
B220, B224, a nd B366 (the last has bee n highUghted and dis- 
cussed also by jPalessandro etal.|2012| because of its anoma- 
lously red {FUV - V)o colour); 

- three metal-rich clusters ([Fe/H]> -1.0) which appear to 
have too blue an HB morphology for their metallicity (B045, 
Gl and B158). In particular, Gl is a confirmed case of 
a metal-rich cluster with both red and blue HB stars (see 
|Meylanetal.|2001| i. 



This feature suggests a similarity with NGC 6441 and 
NGC 6388, the metal-rich Galactic GCs displaying a strong red 
HB clump as well as a significant blue HB tail, that it is currently 
interpreted as due to a large spread in He abundance ( |Caloi &| 
D'Antona||2007| |Yoon et al.|[2U08l ). Also several bright GCs in 



M87 are interpreted as hosting a He abundance spread (see Rey 
|et al.| 2007 their Fig. 9, in particular). Therefore, similarly to 
what is beUeved to happen in some (likely most) Galactic GCs 
(Gratton et al. 2010; Gratton, Carretta, Bragaglia 2012, and ref- 
erences therein), also these M3 1 GCs might host (multiple) pop- 
ulations of He-enriched stars, which evolve as Blue HB stars 
in spite of their high metallicity, as already proposed by, e.g. 



D'Antona et al.,(2005) . The occurrence of the same phenomenon 



in some massive and metal-rich M31 clusters has already been 



suggested with independent observational evidence by Rey et 
|al. (2007); Colucci et al. (2009); Dalessandro et al. | 2012| ). The 
pretty blue {FUV- V)() colour of many Sample B clusters, possi- 
bly indicative of the presence of unseen BHB/EHB populations, 
has been already discussed in Sect. 2.2. 1 above. We stress again 
that Sample B cluster is mainly composed by luminous/massive 
(18 on 25 have My < -9.0), a range of luminosity where only 
a few Galactic GCs (and very peculiar ones, like M54, ll) Cen, 
NGC2419, etc.) are found. 

It may be interesting to explore the available photometric 
and spectroscopic databases in search for further hints on the 
possible existence of multiple stellar populations in M31 GCs. 
Such an analysis is in progress, and the results will be reported 
in a dedicated paper (Fusi Pecci et al., in preparation). 

4. Summary and discussion 

We performed the first thorough and homogeneous comparison 
of the HB morphology between M31 and Milky Way globular 
clusters, based on HB star counts on cluster CMDs. We used 
a simplified version of the Mironov's index (SMI) that has been 
demonstrated to correlate very well with the popular HBR index. 
Strict requirements on the quality of the CMD forced us to limit 
the core of our analysis to a selected sample of 23 M31 GCs, 
mostly located in the outer halo of M31 (Sample A). 

We find that eleven of these clusters lie on a significant dif- 
ferent locus in the metallicity vs HB morphology plane with re- 
spect to their Galactic counterparts located at any distance from 
the MW center. Having an unusually red HB morphology for 
their metallicity, we refer to them as ARSA clusters (see Sect. 
|3]). The possibility that such a difference arose from systematic 
errors in some of the involved parameters has been considered 
and dismissed as very unlikely. 

The most straightforward interpretations of the observed dif- 
ference are in terms of age or He abundance differences: 

1. the ARSA clusters are, on average, » 1 Gyr younger than 
the MW Outer Halo GCs, and ^ 2 Gyr younger than the 
MW inner halo GCs; or 

2. they have a different (internal) distribution of He abundances 
(Y) than their Galactic counterparts, i.e. they should have a 
lower fraction of He-enriched stars (see Sect. [3]). 



Hypothesis 1 seems more natural, by analogy with the Milky 
Way ( Dotter et al. 201 1| . In the Galaxy, slightly younger GCs are 
generally found to be associated with halo substructures and/or 
relics of disrupting dwarf galaxies ( Marin -Franch et aL]|2009| ). 
Also most ARSA clusters can be associated to similar substruc- 
tures in the M31 halo ( [Perina et a l. 2009; Mackey et al.|2010| ), 
and hence younger ages would fit as a general characteristics of 
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recently accreted clusters from dwarf satellite^ The fact that 
Sample A clusters are even younger than the MW OH clus- 
ters may be connected to the more complex and extended re- 
cent accretion activity that took place in Andromeda (w.rt. the 
MW, Ibata et al. 2007| . In the process of hierarchical assem- 
bly of the MW and M31, the most massive (and most metal 
rich) sub-units were likely the first to merge at the bottom of 
the potential well, the strong compression shocks likely lead- 
ing to a large amount of star (and cluster) formation on a short 
timescale. At later times less massive and more metal-poor frag- 
ments (dwarf galaxies) were accreted: these were able to form 
stars over longer timescales (Mateo 1998! Tolstoy, Hill & Tosi 



2009] !. It can be conceived that also the epoch of GC formation 



was consequently postponed in these galaxies. In this context it 
is very interesting to note that, in the [Fe/H] vs. HB morphology 
plane, dwarf galaxy satellites of the MW and M3 1 lie virtually 
in the same locus as ARSA clusters ( Yang & Sarajedini[20l2] ). 

Perhaps the simplest way to envisage hypothesis 2 is to con- 
sider the possibility that the ARSA clusters do not host multi- 
ple populations, as most MW GCs do, or had a less pronounced 
early-enrichment history, hence do not present the sizable He 
abundance spread which acts as an internal second parameter, 
by pushing a fraction of cluster stars to the blue side of the HB. 
It has been suggested that clusters can host multiple populations 
(and consequently He enrichment) only above a certain mass 
threshold (a few 10^ Mq, corresponding to My » -5.5 accord- 
ing to Gratton, Carretta, Bragaglia 2012 and references therein). 
Only one Sample A cluster is below this threshold, while ^ 80% 
of the sample have My < -7.0 and a; 65% are brighter than 
the typical average of GC distributions {My - -7.5, fBrodie & 



Strader, 2006 ), hence they appear as massive as Galactic clus- 
ters where light element (and thus presumably He) enrichment 
has been observed. Moreover, many ARSA clusters present bi- 
modal HB morphology, possibly suggesting that some intrinsic 
difference among cluster stars is indeed there (but see |Fusi Pecci] 



et al. II 1993 for serious warning on simplistic interpretations of 



complex morphologies). The formation environment may also 
have played a role ( Bekki 201 1 1, and the physical conditions that 
led to a later formation epoch may have also influenced the for- 
mation path of these clusters. Therefore, it is also possible that 
both age and He abundance have concurred in producing the ob- 
served diff'erence in HB morphology between outer M31 GCs 
and Milky Way clusters at any distance from the Galactic cen- 
ter It is worth noting that the only Galactic GC in our sample 
that lies on the same isochrone as the ARSA clusters, because 
it has a genuinely redder HB morphology than other clusters of 
the same metallicity (NGC 4590; see Sect.[3]l, is as old as them. 
Therefore, at least in this case, the difference in HB morphol- 
ogy must be due to some other parameter than age (see, e.g. 
|Dalessandroetal.|2012| ). 

Finally, the analysis of sample B clusters provides support 



to the conclusions by Rey et al., (2007 ) and Dalessandro et al. 
2012| l that a population of metal-rich GCs with sizable frac 



tions of Blue HB stars is present in M31, possibly including 
more extreme cases with respect to their Galactic counterparts 
(i.e. NGC 6388, NGC 6441), since they show larger SMI values. 
NGC 6388 and NGC 6441 have a strongly bimodal HB mor- 
phology that is difficult to explain without invoking a significant 
spread in He abundance among cluster stars (Caloi & D'Antona 



' For example, it is interesting to note that the ARSA clusters 
MCGC8 and MCGC9 are both likely associated with the halo filamen- 
tary substructure called Eastern Arc (Stream D, in the nomenclature by 
Mackey et al.|2010^ . 



|2007[ |Yoon et al.]|2008| l. Also, many Sample B clusters have 
{FUV - y)o colours suggesting that the presence of extended 
blue HB tails may be more common in M3 1 clusters than in the 
Milky Way. However, it must be noted that the majority of sam- 
ple B clusters are very bright (72% brighter than My - -9.0 and 
only one cluster with My > -8.0). ICaretta et al.] ( |2010) l found 
that more luminous/massive clusters appear to have wider dis- 
tributions of light elements abundance: this should lead to wider 
distributions of Y and, consequently, to more extended blue tails 
in the HB. Hence, the above result may simply reflect the case 
that very massive clusters, possibly with strong anti-correlation 
patterns and complex HB morphologies including extended blue 
tails (and blue hooks), like co Cen or M54, are much more numer- 
ous in the Andromeda galaxy than in the Milky Wa}{^ Whatever 
the physical mechanism that led to such UV-enhanced HB mor- 
phology, its better knowledge in the M3 1 GC environment could 
certainly add important clues to a better understanding of the 
early evolution of the whole cluster system around our nearby 
galaxy companion, in spite of all the observing problems asso- 
ciated with its large distance. Furthermore, on a larger galactic 
mass scale, this scenario could also shed light to the well recog- 
nized UV-upturn phenomenon, that in most cases seems to con- 



strain the evolution of elliptical galaxies ( Yi et al. 
|& Gonzalez-Lopezlira|2 008 ; Buz zoni et al.|2012 



2011 Buzzoni 
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S. Perina et al.: The horizontal branch morphology of M31 globular clusters. 
Table 1. M3 1 clusters parameters 



Name 


SMI 


[Fe/H] 


E(B-V) 




Rgc 


Sample 










[kpc] 


[kpc] 




BOOS 


0.00±0.10 


-1.00 


0.07 


5.6 


12.1 


A 


B292 


0.91±0.25 


-1.90 


0.15 


16.7 


33.4 


A 


B298 


0.76±0.22 


-1.80 


0.09 


13.9 


28.8 


A 


B336 


0.55±0.26 


-1.90 


0.10 


12.7 


31.5 


A 


B337 


0.00±0.10 


-1.30 


0.06 


13.4 


15.0 


A 


B350 


0.83±0.31 


-1.80 


0.11 


11.4 


15.7 


A 


B407 


0.00±0.10 


-0.40 


0.10 


19.4 


30.9 


A 


B514 


0.61±0.12 


-1.91 


0.09 


54.0 


58.8 


A 


B531 


0.00±0.10 


-0.40 


0.14 


16.8 


35.2 


A 


MCGCl 


0.78±0.22 


-2.15 


0.12 


45.3 


46.9 


A 


MCGC2 


0.77±0.26 


-1.90 


0.10 


32.6 


52.0 


A 


MCGC3 


0.91±0.24 


-1.90 


0.10 


31.1 


39.1 


A 


MCGC5 


0.16±0.05 


-1.90 


0.11 


76.9 


76.9 


A 


MCGC7 


0.00±0.10 


-0.70 


0.06 


17.7 


75.7 


A 


MCGC8 


0.14±0.08 


-1.53 


0.09 


36.2 


54.4 


A 






1 do 


f) 1 f, 

U. 1 u 


Jo.yJ 


fiR 7 
uo.z, 


A 


MCGCl 


0.66±0.18 


-1.90 


0.09 


98.1 


104.0 


A 


MCECl 


0.71±0.14 


-1.91 


0.10 


13.0 


25.2 


A 


MCEC2 


0.11±0.05 


-1.75 


0.13 


35.9 


48.9 


A 


MCEC3 


0.84±0.18 


-1.91 


0.09 


13.7 


13.7 


A 


MCEC4 


0.47±0.18 


-1.78 


0.11 


58.7 


65.4 


A 


B384 


0.00±0.10 


-0.50 


0.04 


16.0 


21.6 


A 


B468 


0.00±0.10 


-0.70 


0.06 


19.7 


22.5 


A 


B006 


0.00±0.10 


-0.55 


0.08 


6.3 


51.4 


B 


BOlO 


0.82±0.26 


-1.80 


0.16 


5.6 


41.5 


B 


B012 


0.94±0.29 


-1.80 


0.11 


5.6 


38.2 


B 


B023 


0.00±0.10 


-0.90 


0.28 


4.3 


54.5 


B 


B027 


0.26±0.15 


-1.66 


0.18 


5.9 


36.6 


B 


B045 


0.39±0.26 


-0.90 


0.16 


4.8 


47.5 


B 


B058 


0.50±0.33 


-1.40 


0.11 


6.8 


25.2 


B 


B088 


0.80±0.21 


-1.90 


0.38 


3.7 


8.0 


B 


B158 


0.61±0.39 


-0.90 


0.09 


2.3 


2.3 


B 


B220 


0.39±0.20 


-1.70 


0.06 


5.1 


22.0 


B 


B224 


0.43±0.19 


-1.80 


0.07 


5.1 


11.8 


B 


B225 


0.00±0.10 


-0.50 


0.05 


4.6 


47.5 


B 


B233 


0.51±0.32 


-1.53 


0.10 


7.9 


8.7 


B 


B240 


0.63+0.37 


-1.66 


0.14 


7.1 


28.5 


B 


B293 


0.93±0.21 


-1.70 


0.12 


16.8 


18.3 


B 


B311 


0.98±0.33 


-1.75 


0.25 


12.8 


39.8 


B 


B338 


0.68±0.20 


-1.20 


0.04 


10.0 


35.8 


B 


B343 


0.86±0.26 


-1.50 


0.10 


14.4 


71.6 


B 


B358 


0.84±0.16 


-1.91 


0.05 


19.4 


76.6 


B 


B366 


0.48±0.11 


-1.80 


0.09 


11.8 


41.6 


B 


B379 


0.00±0.10 


-0.50 


0.13 


11.1 


58.6 


B 


B386 


0.32±0.19 


-1.10 


0.04 


13.8 


19.9 


B 


B405 


0.88±0.20 


-1.55 


0.08 


17.9 


68.7 


B 


GOOl 


0.43±0.09 


-0.90 


0.04 


33.9 


61.9 


B 


B255D 


0.00±0.10 


-0.70 


0.14 


12.3 


18.5 


B 



Notes. The reported uncertainties on SMI are formal propagated Poisson errors on star counts. We arbitrarily assigned an uncertainty of 0. 1 in SMI 
to clusters having only red HB stars. Note that SMI values can be (roughly) converted into HBR, for comparison with other samples, by inverting 
Eq.fT] i.e HBR = IS Ml - 1.0. Distance moduli, [Fe/H], and E(B-V) are taken from Paper-I. Rp is obtained adopting the mean distance modulus 
fromPaper-I for all the clusters, (m - M)o = 24.42. The three-dimensional distance from the center of the galaxy, Rac, is obtained adopting the 
individual distance moduli from Paper-I and (m - M)o = 24.42 for the center of M3 1 . 
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